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Abstract 

Differential inclusive-jet cross sections have been measured for different jet 
radii in neutral current deep inelastic ep scattering for boson virtualities Q 2 > 
125 GeV 2 with the ZEUS detector at HERA using an integrated luminosity of 
81.7 pb _1 . Jets were identified in the Breit frame using the kr cluster algo- 
rithm in the longitudinally inclusive mode for different values of the jet radius 
R. Differential cross sections are presented as functions of Q 2 and the jet trans- 
verse energy, E^ B . The dependence on R of the inclusive-jet cross section has 
been measured for Q 2 > 125 and 500 GeV 2 and found to be linear with R in 
the range studied. Next-to-leading-order QCD calculations give a good descrip- 
tion of the measurements for 0.5<i?<l. A value of a s (Mz) has been extracted 
from the measurements of the inclusive-jet cross-section da/dQ 2 with R = 1 for 
Q 2 > 500 GeV 2 : a s {M z ) = 0.1207 ± 0.0014 (stat.) t°;°°l (exp.) i°;°° 22 (th.). 
The variation of a s with E^ B is in good agreement with the running of a s as 
predicted by QCD. 
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1 Introduction 



The study of jet production in ep collisions at HERA has been well established as a 
testing ground of perturbative QCD (pQCD) providing precise determinations of the 
strong coupling constant, a s , and its scale dependence. The jet observables used to 
test pQCD included dijet [1-4], inclusive-jet [2,4-6] and multijet [7,8] cross sections in 
neutral current (NC) deep inelastic ep scattering (DIS), dijet [9-13], inclusive-jet [14,15] 
and multijet [16] cross sections in photoproduction and the internal structure of jets in 
NC [17-19] and charged current [20] DIS. 

These studies demonstrated that the cluster algorithm [21] in the longitudinally in- 
variant inclusive mode [22] is at present the method to reconstruct jets in ep collisions 
for which the smallest uncertainties are achieved. Previous analyses were done with jet 
radius R = 1, where R is the maximum distance in the pseudorapidity (77) - azimuth 
(4>) plane for particle recombination. The study of the predicted jet cross sections us- 
ing different jet radii allows the identification of the values of R for which the theory is 
most reliable. Furthermore, smaller values of the jet radius R are of particular interest 
for the identification of heavy particles decaying into jets; in these decays, the final-state 
jets may emerge close in phase space and need to be identified separately for a faithful 
reconstruction of the properties of the parent particle [23]. Neutral current DIS provides 
a well understood environment in which to study the dependence of jet production on 
the jet radius and to confront the data with precise next-to-leading-order (NLO) QCD 
calculations in a hadron- induced reaction. 

The hadronic final state in NC DIS may consist of jets of high transverse energy, E^*, 
produced in the hard-scattering process. In NC DIS, the Breit frame [24] is preferred to 
conduct the jet search, since jet production is directly sensitive to hard QCD processes; 
in the Born process (eg — > eg), the virtual boson (V*, with V* = j,Z) is absorbed by 
the struck quark, which is back-scattered with zero transverse momentum with respect 
to the V* direction. At leading order (LO) in a s , the boson-gluon-fusion (V*g — > qq) 
and QCD-Compton (V*q — > qg) processes give rise to two hard jets with opposite trans- 
verse momenta. Inclusive-jet production at high allows more stringent tests of the 
pQCD calculations than dijet production due to the reduced theoretical uncertainties. 
Restrictions on the topology of dijet events are necessary to avoid infrared-sensitive re- 
gions where the NLO QCD programs are not reliable, whereas inclusive-jet cross sections 
are infrared insensitive. Therefore, such measurements allow tests of pQCD in the widest 
phase-space region for jet production. In particular, previous measurements of inclusive- 
jet cross sections in NC DIS at high Q 2 [5], where Q 2 is the negative of the square of 
the four-momentum transfer, provided the most precise determination of a s at HERA to 
date. 
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This letter presents new measurements of differential inclusive-jet cross sections as a 
function of the jet transverse energy in the Breit frame, E^ B , and Q 2 for different values 
of R. For R — 1, this analysis is based on the same inclusive-jet data sample presented in a 
recent publication [4]. The results have been compared with NLO QCD calculations using 
recent parameterisations of the parton distribution functions (PDFs) of the proton [25-27]. 
In addition, an updated determination of a s and of its scale dependence has been obtained 
using a data sample which corresponds to more than a twofold increase in luminosity with 
respect to the previous analysis [5]. 



2 Experimental set-up 

A detailed description of the ZEUS detector can be found elsewhere [28,29]. A brief 
outline of the components that are most relevant for this analysis is given below. 

Charged particles are tracked in the central tracking detector (CTD) [30], which operates 
in a magnetic field of 1.43 T provided by a thin superconducting solenoid. The CTD 
consists of 72 cylindrical drift-chamber layers, organised in nine superlayers covering the 
polar-angle 1 region 15° < 9 < 164°. The transverse-momentum resolution for full-length 
tracks can be parameterised as u(pt)/pt = 0.0058j>r © 0.0065 © 0.0014/pr, with p T in 
GeV. The tracking system was used to measure the interaction vertex with a typical 
resolution along (transverse to) the beam direction of 0.4 (0.1) cm and to cross-check the 
energy scale of the calorimeter. 

The high-resolution uranium-scintillator calorimeter (CAL) [31] covers 99.7% of the total 
solid angle and consists of three parts: the forward (FCAL), the barrel (BCAL) and 
the rear (RCAL) calorimeters. Each part is subdivided transversely into towers and 
longitudinally into one electromagnetic section and either one (in RCAL) or two (in BCAL 
and FCAL) hadronic sections. The smallest subdivision of the calorimeter is called a 
cell. Under test-beam conditions, the CAL single-particle relative energy resolutions were 
a(E)/E = 0.18/ VE for electrons and a(E)/E = 0.35/ VE for hadrons, with E in GeV. 

The luminosity was measured from the rate of the bremsstrahlung process ep — > e^p. The 
resulting small-angle energetic photons were measured by the luminosity monitor [32], a 
lead-scintillator calorimeter placed in the HERA tunnel at Z = —107 m. 



The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. 
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3 Data selection and jet search 



The data were collected during the running period 1998-2000, when HERA operated with 
protons of energy E p = 920 GeV and electrons or positrons 2 of energy E e = 27.5 GeV, and 
correspond to an integrated luminosity of 81.7 ± 1.8 pb _1 , of which 16.7 pb" 1 (65.0 ptr 1 ) 
was for e~p (e + p) collisions. 

Neutral current DIS events were selected offline using the same criteria as reported in a 
recent publication [4]. The main steps are briefly listed below. 

The scattered-electron candidate was identified from the pattern of energy deposits in the 
CAL [33]. The energy (E' e ) and polar angle (9 e ) of the electron candidate were determined 
from the CAL measurements. The Q 2 variable was reconstructed using the double angle 
method (Qda) [34]- The angle jh, which corresponds to the angle of the scattered quark 
in the quark-parton model, was reconstructed using the hadronic final state [34]. 

The main requirements imposed on the data sample were: an electron candidate with 
E' e > 10 GeV; a vertex position along the beam axis in the range \Z\ < 34 cm; 38 < 
(E — Pz) < 65 GeV, where E is the total energy as measured by the CAL, E = 
and Pz is the Z-component of the vector P = J^-Eir;; in both cases the sum runs over 
all CAL cells, Ei is the energy of the CAL cell % and rj is a unit vector along the line 
joining the reconstructed vertex and the geometric centre of the cell i; Q 2 > 125 GeV 2 ; 
and | cos7h| < 0.65. After all these requirements, contamination from non-ep interactions 
and other physics processes was negligible. 

The kx cluster algorithm was used in the longitudinally invariant inclusive mode to recon- 
struct jets in the hadronic final state both in data and in Monte Carlo (MC) simulated 
events (see Section 4). In the data, the algorithm was applied to the energy deposits 
measured in the CAL cells after excluding those associated with the scattered-electron 
candidate. In the following discussion, E % TB denotes the transverse energy, rj B the pseu- 
dorapidity and <p B the azimuthal angle of object % in the Breit frame. For each pair of 
objects, where the initial objects are the energy deposits in the CAL cells, the quantity 

da = min(^ jB , E 3 TB ) 2 ■ [( Vb - Vb ) 2 + {<& - ^ B ) 2 ]/R 2 

was calculated. For each individual object, the quantity di = (E^ B ) 2 was also calculated. 
If, of all the values {dij,di}, d^i was the smallest, then objects k and / were combined 
into a single new object. If, however, dk was the smallest, then object k was considered a 
jet and was removed from the sample. The procedure was repeated until all objects were 

2 Here and in the following, the term "electron" denotes generically both the electron (e~) and the 
positron (e + ). 
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assigned to jets. The jet search was performed with different jet radii (R = 0.5, 0.7 and 
1) and the jet variables were defined according to the Snowmass convention [35]. 

After reconstructing the jet variables in the Breit frame, the massless four-momenta were 
boosted into the laboratory frame, where the transverse energy (E 3 ^ hAB ) and the pseudo- 
rapidity (?7lab) °f eacn were calculated. Energy corrections [5,11,18] were then applied 
to the jets in the laboratory frame and propagated back into E^ B . The jet variables in 
the laboratory frame were also used to apply additional cuts on the selected sample, as 
explained in a recent publication [4]: events were removed from the sample if any of the 
jets was in the backward region of the detector (t/lab < ~ 2) and jets were not included 
in the final sample if -E^lab < 2.5 GeV. 

Only events with at least one jet in the pseudorapidity range —2 < 77^ < 1.5 were kept 
for further analysis. The final data samples with at least one jet satisfying E 3 ^ > 8 GeV 
contained 19908 events for R = 1, 16231 for R = 0.7 and 12935 for R = 0.5. 



4 Monte Carlo simulation 

Samples of events were generated to determine the response of the detector to jets of 
hadrons and the correction factors necessary to obtain the hadron-level jet cross sections. 
The hadron level is defined by those hadrons with lifetime r > 10 ps. The generated events 
were passed through the Geant 3.13-based [36] ZEUS detector- and trigger-simulation 
programs [29]. They were reconstructed and analysed by the same program chain as the 
data. 

Neutral current DIS events including electroweak radiative effects were simulated using 
the Heracles 4.6.1 [37] program with the Djangoh 1.1 [38] interface to the QCD 
programs. The QCD cascade is simulated using the colour-dipole model (CDM) [39] in- 
cluding the leading-order (LO) QCD diagrams as implemented in Ariadne 4.08 [40] and, 
alternatively, with the MEPS model of Lepto 6.5 [41]. The CTEQ5D [42] parameterisa- 
tions of the proton PDFs were used for these simulations. Fragmentation into hadrons is 
performed using the Lund string model [43] as implemented in Jetset [44,45]. 

The jet search was performed on the MC events using the energy measured in the CAL 
cells in the same way as for the data. The same jet algorithm was also applied to the 
final-state particles (hadron level) and to the partons available after the parton shower 
(parton level). The MC programs were also used to correct the measured cross sections 
for QED radiative effects and the running of a cm . 



4 



5 QCD calculations 



The measurements were compared with LO (0(a s )) and NLO QCD (0(a 2 )) calculations 
obtained using the program Disent [46]. The calculations were performed in the MS 
renormalisation and factorisation schemes using a generalised version [46] of the subtrac- 
tion method [47]. The number of flavours was set to five and the renormalisation (/xr) 
and factorisation (jif) scales were chosen to be [Mr = E^ B and /ip = Q, respectively. The 

(5) 

strong coupling constant was calculated at two loops with A^± = 226 MeV, corresponding 
to a s (Mz) = 0.118. The calculations were performed using the ZEUS-S [25] parameter- 
isations of the proton PDFs. In Disent, the value of a cm was fixed to 1/137. The 
cluster algorithm was also applied to the partons in the events generated by Disent in 
order to compute the jet cross-section predictions. The calculations were performed for 
the same values of R as for the data. In addition, predictions were obtained for R = 0.3 
and 1.2 to determine the range of R in which the theory is most reliable. 

Since the measurements refer to jets of hadrons, whereas the NLO QCD calculations refer 
to jets of partons, the predictions were corrected to the hadron level using the MC models. 
The multiplicative correction factor (Chad) was defined as the ratio of the cross section 
for jets of hadrons over that for jets of partons, estimated by using the MC programs 
described in Section 4. The mean of the ratios obtained with Ariadne and Lepto- 
MEPS was taken as the value of Chad- The value of Chad differs from unity by less than 
5%, 15% and 25% for R = 1, 0.7 and 0.5, respectively, in the region Q 2 > 500 GeV 2 . For 
R = 1.2, Chad differs from unity by less than 1% and for R = 0.3, Chad differs by 40%. 

Disent does not include the contribution from Z° exchange; MC simulated events with 
and without Z° exchange were used to include this effect in the NLO QCD predictions. In 
the following, NLO QCD calculations will refer to the fully corrected predictions, except 
where otherwise stated. 

5.1 Theoretical uncertainties 

Several sources of uncertainty in the theoretical predictions were considered: 

• the uncertainty on the NLO QCD calculations due to terms beyond NLO, estimated 
by varying /jr between E 3 ^ t B /2 and 2£j? t B , was below ±7% at low Q 2 and low E } ^ B 
and decreased to less than ±5% for Q 2 > 250 GeV 2 for R — 1. For smaller radii, 
the estimated uncertainty is smaller (higher) at low (high) Q 2 than for R — 1. For 
R = 1.2, this uncertainty increases up to ±10% for Q 2 rs 500 GeV 2 ; 

• the uncertainty on the NLO QCD calculations due to those on the proton PDFs was 
estimated by repeating the calculations using 22 additional sets from the ZEUS-S anal- 
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ysis, which takes into account the statistical and correlated systematic experimental 
uncertainties of each data set used in the determination of the proton PDFs. The 
resulting uncertainty in the cross sections was below ±3%, except in the high-.Ej? B 
region where it reached ±4.4%; 

• the uncertainty on the NLO QCD calculations due to that on a s [Mz) was estimated by 
repeating the calculations using two additional sets of proton PDFs, for which different 
values of a s (Mz) were assumed in the fits. The difference between the calculations 
using these various sets was scaled by a factor such as to reflect the uncertainty on 
the current world average of a s [48]. The resulting uncertainty in the cross sections 
was below ±2%; 

• the uncertainty from the modelling of the QCD cascade was assumed to be half the 
difference between the hadronisation corrections obtained using the Ariadne and 
Lepto-MEPS models. The resulting uncertainty on the cross sections was less than 
1.4% for R = 1 and increased up to ~ 4% for R = 0.5; 

• the uncertainty of the calculations in the value of \ip was estimated by repeating the 
calculations with \ip = Q/2 and 2Q. The variation of the calculations was negligible. 

The total theoretical uncertainty was obtained by adding in quadrature the individual 
uncertainties listed above. 

It is concluded that NLO QCD provides predictions with comparable precision in the range 
R = 0.5 — 1. For larger values of R, e.g. R = 1.2, it was estimated that the uncertainty 
on the NLO QCD calculations due to terms beyond NLO increases up to about 10% for 
high Q 2 values. On the other hand, the hadronisation correction estimated for the cross 
sections with smaller radii, e.g. R = 0.3, increases up to about 40%. Therefore, only 
measurements for the range R = 0.5 — 1 are presented in Section 7. 

6 Acceptance corrections 

The E^ B and Q 2 distributions in the data were corrected for detector effects using bin-by- 
bin correction factors determined with the MC samples. These correction factors took into 
account the efficiency of the trigger, the selection criteria and the purity and efficiency of 
the jet reconstruction. For this approach to be valid, the uncorrected distributions of the 
data must be well described by the MC simulations at the detector level. This condition 
was satisfied by both the Ariadne and Lepto-MEPS MC samples. The average between 
the acceptance-correction values obtained with Ariadne and Lepto-MEPS was used to 
correct the data to the hadron level. The deviations in the results obtained by using either 
Ariadne or Lepto-MEPS to correct the data from their average were taken to represent 



6 



systematic uncertainties of the effect of the QCD-cascade model in the corrections (see 
Section 6.1). The acceptance-correction factors differed from unity by typically less than 
10%. 



6.1 Experimental uncertainties 

The following sources of systematic uncertainty were considered for the measured cross 
sections: 

• the uncertainty in the absolute energy scale of the jets was estimated to be ±1% 
for E^ LAB > 10 GeV and ±3% for lower E 3 ^ LAB values [10,11,49]. The resulting 
uncertainty was about ±5%; 

• the differences in the results obtained by using either Ariadne or Lepto-MEPS to 
correct the data for detector effects were taken to represent systematic uncertainties. 
The resulting uncertainty was typically below ±3%; 

• the uncertainty due to the selection cuts was estimated by varying the values of the 
cuts within the resolution of each variable; the effect on the cross sections was typically 
below ±3%; 

• the uncertainty on the reconstruction of the boost to the Breit frame was estimated 
by using the direction of the track associated to the scattered electron instead of that 
derived from its impact position in the CAL. The effect was typically below ±1%; 

• the uncertainty in the absolute energy scale of the electron candidate was estimated 
to be ±1% [50]. The resulting uncertainty was below ±1%; 

• the uncertainty in the cross sections due to that in the simulation of the trigger was 
negligible. 

The systematic uncertainties not associated with the absolute energy scale of the jets were 
added in quadrature to the statistical uncertainties and are shown in the figures as error 
bars. The uncertainty due to the absolute energy scale of the jets is shown separately as a 
shaded band in each figure, due to the large bin-to-bin correlation. In addition, there was 
an overall normalisation uncertainty of 2.2% from the luminosity determination, which is 
not included in the figures. 
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7 Results 



7.1 Inclusive-jet differential cross sections for different jet radii 

The inclusive-jet differential cross sections were measured in the kinematic region Q 2 > 
125 GeV 2 and | cos7/j| < 0.65. These cross sections include every jet of hadrons in the 
event with E^ B > 8 GeV and —2 < r/ 3 ^ < 1.5 and were corrected for detector and QED 
radiative effects and the running of a em . 

The measurements of the inclusive-jet differential cross sections as functions of E 3 ^ B and 
Q 2 are presented in Figs, la and 2a for jet radii R = l, 0.7 and 0.5. In these figures, each 
data point is plotted at the abscissa at which the NLO QCD differential cross section 
was equal to its bin-averaged value. The measured da/dE 3 ^ (da/dQ 2 ) exhibits a steep 
fall-off over three (five) orders of magnitude for the jet radii considered in the E^ B (Q 2 ) 
measured range. 

The NLO QCD predictions with jin = E 3 ^ are compared to the measurements in Figs, la 
and 2a. The fractional difference of the measured differential cross sections to the NLO 
QCD calculations is shown in Figs, lb and 2b. The calculations reproduce the measured 
differential cross sections well for the jet radii considered, with similar precision. To study 
the scale dependence, NLO QCD calculations using /ir = Q were also compared to the 
data (not shown); they also provide a good description of the data. 

7.2 Dependence of the inclusive-jet cross section on the jet ra- 
dius 

Measurements of the inclusive-jet cross section have been performed for > 8 GeV 
and —2 < rf B < 1.5 in the kinematic range given by | cos7h| < 0.65 integrated above 
Qmm = 125 and 500 GeV 2 for different jet radii. The measured cross section, <Xj ets , as 
a function of R is presented in Figs. 3a and 3b. The measured cross sections increase 
linearly with R in the range between 0.5 and 1. The increase of CTj Cts as R increases can 
be understood as the result of more transverse energy being gathered in a jet so that a 
larger number of jets has E 3 ^ exceeding the threshold of 8 GeV. 

The predictions of LO and NLO QCD for Oj Cts at the parton level, with no corrections for 
hadronisation or Z° exchange, are shown in the inset of Fig. 3. The LO predictions do 
not depend on R since there is only one parton per jet. The NLO calculations give the 
lowest-order contribution to the R dependence of the inclusive-jet cross section. The NLO 
QCD calculations, corrected to include hadronisation and Z° effects, are also shown in 
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the inset of Fig. 3. It is observed that, in addition to the effects due to parton radiation, 
the hadronisation corrections further modify the shape of the prediction. 

The NLO QCD calculations are compared to the data in Fig. 3. They give a good 
description of the R dependence of the data within the jet-radius range considered. The 
uncertainty of the NLO calculation is also shown. The total theoretical uncertainty of the 
cross section integrated above Q^in = 125 GeV 2 (500 GeV 2 ) changes from 5.6% (3.2%) 
for R = 1 to 4.2% (7.1%) for R = 0.5. This variation is due to a change of the uncertainty 
arising from higher-order QCD corrections and hadronisation. The uncertainties arising 
from a s [Mz) and the proton PDFs do not change significantly with R. 

7.3 Determination of a. s (Mz) 

The measured differential cross sections presented in Section 7.1 were used to determine a 
value of a s {Mz) using the method presented previously [5]. The NLO QCD calculations 
were performed using the program Disent with five sets of ZEUS-S proton PDFs which 
were determined from global fits assuming different values of a s (Mz), namely a s (Mz) = 
0.115, 0.117, 0.119, 0.121 and 0.123. The value of a s [Mz) used in each calculation was 
that associated with the corresponding set of PDFs. The a s [Mz) dependence of the 
predicted cross sections in each bin % of A (A = Q 2 , E 3 ^) was parameterised according to 

[da/dA(a a (M z ))]i = C\a s {M z ) + C\c? s {Mz), 

where C\ and C\ were determined from a x 2 fit to the NLO QCD calculations. The value 
of a s (Mz) was determined by a \ 2 fit to the measured da/dA values for several regions of 
the variable A. The values of a s (Mz) obtained from the various differential cross sections 
and jet radii were found to be consistent. The result obtained using the measured da/dQ 2 
for Q 2 > 500 GeV 2 with R = 1 yields the smallest uncertainty and, therefore, represents 
the most precise determination from this analysis. In the region Q 2 > 500 GeV 2 , the 
experimental systematic uncertainties on the cross sections are smaller than at lower Q 2 
and the theoretical uncertainties due to the proton PDFs and to terms beyond NLO are 
minimised. 

The uncertainties on the extracted values of a s (Mz) due to the experimental systematic 
uncertainties were evaluated by repeating the analysis for each systematic check presented 
in Section 6.1. The overall normalisation uncertainty from the luminosity determination 
was also considered. The largest contribution to the experimental uncertainty comes 
from the jet energy scale and amounts to ±2% on a s (Mz)- The theoretical uncertainties 
were evaluated as described in Section 5.1. The largest contribution was the theoretical 
uncertainty arising from terms beyond NLO, which was estimated by using the method 
proposed by Jones et al. [51], and amounted to ±1.5%. The uncertainty due to the proton 
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PDFs was ±0.7%. The uncertainty arising from the hadronisation effects amounted to 
±0.8%. 

As a cross-check, a s (Mz) wa s determined by using NLO QCD calculations based on the 
MRST2001 [26] and CTEQ6 [27] sets of proton PDFs. The values obtained are consistent 
within 1% with those based on ZEUS-S. The uncertainty arising from the proton PDFs was 
estimated to be ±0.7% (±1.6%) using the results of the MRST2001 (CTEQ6) analysis. 

The value of a s {Mz) obtained from the measured da/dQ 2 for Q 2 > 500 GeV 2 with R = 1 
is 

a s (M z ) = 0.1207 ±0.0014 (stat.) t°;°ol (exp.) ±H°23 (th-)- 

This value of a s (Mz) is consistent with the current world average of 0.1189 ± 0.0010 [48] 
as well as with the HERA average of 0.1186 ± 0.0051 [52]. It has a precision comparable 
to the values obtained from e + e~ interactions [48]. 

7.4 Energy-scale dependence of a s 

The QCD prediction for the energy-scale dependence of the strong coupling constant was 
tested by determining a s from the measured da/dE^ with R — 1 at different E^ B 
values. The method employed was the same as described above, but parameterising the 
a s dependence of da/dE^ in terms of a s ((i?y 1 B )) instead of a s (Mz), where (-E^b) is the 
average of the data in each bin. The extracted values of a s are shown in Fig. 4. 
The results are in good agreement with the predicted running of the strong coupling 
constant [53] calculated at two loops [54] over a large range in £y Ct B . 



8 Summary 

Measurements of differential cross sections for inclusive-jet production in neutral current 
deep inelastic ep scattering at a centre-of-mass energy of 318 GeV have been presented. 
The cross sections refer to jets of hadrons identified in the Breit frame with the kr cluster 
algorithm in the longitudinally invariant inclusive mode. The cross sections are given in 
the kinematic region of Q 2 > 125 GeV 2 and | cos7^| < 0.65. 

The dependence of the inclusive-jet cross sections on the jet-radius R has been studied. 
It has been determined that NLO QCD provides predictions with comparable precision 
in the range R = 0.5 — 1. Measurements of inclusive-jet differential cross sections for 
this jet-radius range have been presented. The NLO QCD calculations provide a good 
description of the measured inclusive-jet differential cross sections da/dE 3 ^ B and da/dQ 2 
for R = 1, 0.7 and 0.5. It is observed that the measured inclusive-jet cross section 
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integrated above Q^ in = 125 and 500 GeV 2 increases linearly with R in the jet-radius 
range studied. 

The measured inclusive-jet differential cross sections have been used to extract a value of 
a s {M z ). A QCD fit of the cross-section da/dQ 2 with R = 1 for Q 2 > 500 GeV 2 yields 
the determination with smallest uncertainty, 

a s (M z ) = 0.1207± 0.0014 (stat.) l^ol (exp.) Ijffi (th.). 

This value is in good agreement with the world and HERA averages. The extracted values 
of a s at different * B are in good agreement with the predicted running of the strong 
coupling constant over a large range in E^ B . 
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Ey* B bin 
(GeV) 



da/dE^ 
(pb/GeV) 



J stat 



-'syst 



a 



QED 



R 



8-10 


63.98 


0.68 


10-14 


29.29 


0.34 


14-18 


11.07 


0.20 


18-25 


3.234 


0.080 


25-35 


0.773 


0.033 


35-100 


0.0312 


0.0027 



+1.20 


+2.84 


-1.39 


-2.56 


+0.37 


+1.25 


-0.52 


-1.31 


+0.15 


+0.61 


-0.20 


-0.50 


+0.036 


+0.156 


-0.045 


-0.167 


+0.016 


+0.040 


-0.015 


-0.033 


+0.0005 


+0.0015 


-0.0006 


-0.0022 



0.95 
0.96 
0.96 
0.94 
0.95 
1.06 



R = 0.7 



8-10 


50.09 


0.60 


10-14 


23.38 


0.30 


14-18 


8.97 


0.18 


18-25 


2.659 


0.071 


25-35 


0.631 


0.029 


35-100 


0.0237 


0.0022 



+0.93 


+2.16 


-1.19 


-1.94 


+0.30 


+1.11 


-0.38 


-1.06 


+0.17 


+0.47 


-0.19 


-0.41 


+0.031 


+0.120 


-0.038 


-0.126 


+0.024 


+0.028 


-0.025 


-0.026 


+0.0006 


+0.0015 


-0.0007 


-0.0016 



0.95 
0.96 
0.96 
0.95 
0.96 
1.03 



R = 0.5 



8-10 


38.25 


0.52 


10-14 


17.78 


0.26 


14-18 


7.09 


0.15 


18-25 


2.257 


0.063 


25-35 


0.514 


0.025 


35-100 


0.0208 


0.0019 



+0.85 


+1.81 


-1.11 


-1.60 


+0.26 


+0.85 


-0.30 


-0.77 


+0.22 


+0.32 


-0.23 


-0.33 


+0.027 


+0.096 


-0.032 


-0.103 


+0.020 


+0.024 


-0.020 


-0.024 


+0.0006 


+0.0014 


-0.0006 


-0.0011 



0.95 
0.96 
0.95 
0.96 
0.97 
1.04 



had 



0.91 
0.95 
0.96 
0.97 
0.95 
0.95 



0.77 
0.83 
0.88 
0.91 
0.92 
0.93 



0.64 
0.70 
0.77 
0.83 
0.86 
0.90 



Table 1: Inclusive jet cross- sections da/dE 3 ^ for jets of hadrons in the Breit 
frame selected with the longitudinally invariant kT cluster algorithm for different 
values of R (Fig. 1). The statistical, uncorrelated systematic and jet- energy- scale 
(ESj uncertainties are shown separately. The multiplicative corrections applied 
to the data to correct for QED radiative effects, Cqed; and the corrections for 
hadronisation effects to be applied to the parton-level NLO QCD calculations, Chad? 
are shown in the last two columns. 
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Q 2 bin 


da/dQ 2 












(GeV 2 ) 


(pb/GeV 2 ) 


^stat 


^syst 




Cqed 


Chad 






R = 


: 1 








125-250 


1.106 


0.012 


+0.013 
-0.020 


+0.066 
-0.062 


0.97 


0.92 


250-500 


0.3671 


0.0053 


+0.0048 
-0.0078 


+0.0153 
-0.0149 


0.95 


0.94 


500-1000 


0.1037 


0.0020 


+0.0020 
-0.0021 


+0.0033 
-0.0029 


0.95 


0.95 


1000-2000 


0.02439 


0.00072 


+0.00039 
-0.00033 


+0.00059 
-0.00058 


0.94 


0.96 


2000-5000 


0.00396 


0.00017 


+0.00017 
-0.00015 


+0.00008 
-0.00008 


0.94 


0.95 


5000-100000 


0.000036 


0.000003 


+0.000003 
-0.000003 


+0.000001 
-0.000001 


0.98 


0.96 






R = 


0.7 








125-250 


0.855 


0.010 


+0.007 
-0.012 


+0.054 
-0.048 


0.97 


0.79 


250-500 


0.2913 


0.0046 


+0.0053 
-0.0079 


+0.0124 
-0.0119 


0.95 


0.83 


500-1000 


0.0840 


0.0018 


+0.0017 
-0.0017 


+0.0026 
-0.0024 


0.95 


0.86 


1000-2000 


0.02079 


0.00066 


+0.00041 
-0.00043 


+0.00046 
-0.00049 


0.94 


0.88 


2000-5000 


0.00332 


0.00016 


+0.00018 
-0.00016 


+0.00007 
-0.00006 


0.93 


0.88 


5000-100000 


0.000031 


0.000003 


+0.000002 
-0.000002 


+0.000001 
-0.000001 


0.97 


0.90 






R = 


0.5 








125-250 


0.6344 


0.0088 


+0.0058 
-0.0092 


+0.0406 
-0.0357 


0.97 


0.64 


250-500 


0.2246 


0.0040 


+0.0053 
-0.0069 


+0.0097 
-0.0097 


0.95 


0.70 


500-1000 


0.0672 


0.0016 


+0.0020 
-0.0019 


+0.0021 
-0.0020 


0.94 


0.75 


1000-2000 


0.01709 


0.00060 


+0.00049 
-0.00051 


+0.00043 
-0.00042 


0.94 


0.79 


2000-5000 


0.00296 


0.00015 


+0.00016 
-0.00015 


+0.00006 
-0.00006 


0.95 


0.81 


5000-100000 


0.000028 


0.000003 


+0.000002 
-0.000002 


+0.000001 
-0.000000 


0.98 


0.83 



Table 2: Inclusive jet cross-sections da/dQ 2 for jets of hadrons in the Breit frame 
selected with the longitudinally invariant cluster algorithm for different values 
of R (Fig. 2). Other details as in the caption to Table 1. 
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R 


(Tjets 

(Pb) <W 4yst ^ES 


Cqed 


Chad 


Q 2 > 125 GeV 2 


0.5 
0.7 
1.0 


197.8 1.9 tli tU 
255.6 2.1 tU t\{i 
321.5 2.4 tU t\ii 


0.96 
0.96 
0.96 


0.70 
0.82 
0.94 


Q 2 > 500 GeV 2 


0.5 
0.7 
1.0 


62.3 1.1 IJj t\]j 
75.8 1.3 

91 6 14 +L6 + 2 - 6 
yi.u i.i _ 15 _ 24 


0.95 
0.95 
0.95 


0.77 
0.87 
0.95 



Table 3: Inclusive jet cross- sections <7j ets for jets of hadrons in the Breit frame 
selected with the longitudinally invariant kx cluster algorithm for Q 2 > 125 and 
500 GeV 2 (Fig. 3). Other details as in the caption to Table 1. 



(4%> 


ct s 








(GeV) 




Ostat 


S yst 


^theor 


8.9 


0.1907 


+0.0038 
-0.0038 


+0.0194 
-0.0171 


+0.0208 
-0.0192 


11.7 


0.1746 


+0.0028 
-0.0028 


+0.0123 
-0.0126 


+0.0148 
-0.0142 


15.7 


0.1719 


+0.0032 
-0.0031 


+0.0105 
-0.0092 


+0.0107 
-0.0105 


20.7 


0.1519 


+0.0028 
-0.0028 


+0.0061 
-0.0065 


+0.0057 
-0.0057 


28.6 


0.1512 


+0.0037 
-0.0037 


+0.0050 
-0.0045 


+0.0043 
-0.0044 


41.2 


0.1452 


+0.0064 
-0.0063 


+0.0041 
-0.0056 


+0.0036 
-0.0036 



Table 4: The a s values determined from a Q CD fit of the measured da/dE s ^ B with 

R = 1 as a function of E 3 ^ B (Fig. 4)- The statistical, systematic and theoretical 
uncertainties are shown separately. 



17 



ZEUS 




Figure 1: (a) The measured differential cross-section da/dE 3 ^ B for inclusive-jet 

production with —2 < rf B < 1.5 (dots) for different jet radii, in the kinematic range 
given by Q 2 > 125 GeV 2 and |cos7/j| < 0.65. The NLO QCD calculations with 
[ir = -E^b (solid lines), corrected to include hadronisation and Z° effects and using 
the ZEUS-S parameterisations of the proton PDFs, are also shown. Each cross 
section has been multiplied by the scale factor indicated in brackets to aid visibility, 
(b) The fractional differences between the measured da/dE 3 ^ and the NLO QCD 
calculations ( dots ); the hatched bands display the total theoretical uncertainty. The 
inner error bars represent the statistical uncertainty. The outer error bars show 
the statistical and systematic uncertainties, not associated with the uncertainty in 
the absolute energy scale of the jets, added in quadrature. The shaded bands display 
the uncertainty due to the absolute energy scale of the jets. 
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Figure 2: T7ie measured differential cross-section da/dQ 2 for inclusive-jet pro- 
duction with Ej? B > 8 GeV and —2 < < 1.5 (dots) for different jet radii, in the 
kinematic range given by | COS7/J < 0.65. Other details as in the caption to Fig. 1. 
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Figure 3: The measured cross-section o~j Cts as a function of the jet radius for 
inclusive-jet production with E^ B > 8 GeV and —2 < rf^ < 1.5 (dots), in the 
kinematic range given by 10087^1 < 0.65 and (a) Q 2 > 125 GeV 2 and (b) Q 2 > 
500 GeV 2 . The insets show the LO (dot-dashed lines) and NLO (dashed lines) 
QCD calculations. The NLO QCD calculations corrected to include hadronisation 
and Z° effects are shown as solid lines. Other details as in the caption to Fig. 1. 
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Figure 4: The a s values determined from the measured da/dE s ^ B with R = 1 
as a function of E^ B (dots). The dashed line indicates the renormalisation group 
prediction at two loops obtained from the a s (Mz) value determined in this analysis 
and the shaded area represents its uncertainty. The inner error bars represent the 
statistical uncertainties of the data. The outer error bars show the statistical and 
systematic uncertainties added in quadrature. The dotted vertical bars, shifted to 
aid visibility, represent the theoretical uncertainties. 
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